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MULTI- FREQUENCY TONE DETECTOR 
Field Of The Invention 

The present invention relates generally to tone detectors for telephone systems 
and in particular to an improved tone detector with filter characterized by spectral 
5 nulls located to suppress frequencies adjacent the tone to be detected and an improved 
detection algorithm using roll-off characteristics of the filter and background noise. 

Background Of the Invention 

10 Multi-Frequency (MF) tones are used extensively in telecommunications. 

Examples of MF tones are DTMF and R2. Typical prior art tone detectors include a 
pre-filter for suppressing out-of-band frequencies, an analysis filter for determining 
the energy of the input signal at each multi-frequency tone of interest, and a decision 
logic block for interpreting the multi-frequency energy values output from the 

15 analysis filter and in response performing a pass/fail test for determining the presence 
or absence of a tone at each multi-frequency tone of interest. 

Prior art analysis filters have been implemented using standard digital filtering 
techniques. The analysis window for such prior art digital filters must be carefully 
20 chosen. Rectangular windows are often used for simplicity. However, the leakage of 
energies in the side lobes of the rectangular window is such that other MF frequencies 
may be incorrectly detected. Also, fast roll-off characteristics (i.e. high order) are 
usually imposed in order to fail any tone that deviates excessively from its nominal 
value. 

25 

Summary Of The Invention 

According to one aspect of the present invention there is provided a multi- 
frequency tone detector with analysis window (i.e. filter size) chosen such that 
30 spectral nulls are located at adjacent frequencies of interest. According to another 
aspect, the decision logic block uses the roll-off characteristics of the filter in 
conjunction with background noise to determine a threshold pass/fail for any tone that 
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has deviated excessively from its nominal value. The foregoing aspects of the 
invention result in simpler filter design (i.e. reduced order) relative to prior art tone 
detectors, while maintaining high tone detection accuracy. 

5 Brief Description Of The Drawings 

A preferred embodiment of the present invention will now be described more 
fully with reference to the accompanying drawings in which: 

10 Figure 1 is a block diagram of a well known multi-frequency tone detector; 

Figure 2 shows the frequency response of a tone detector filter with 
rectangular window; 

15 Figure 3 shows the frequency response of a filter with rectangular window 

centered at 900 Hz, showing the location of spectral nulls at equally spaced 
frequencies adjacent 900 Hz, in accordance with the present invention; 



Figure 5 shows the frequency response of a filter with rectangular window of 
size N— 195, centered at 852 Hz, showing moderate suppression of adjacent DTMF 
signal frequencies, in accordance with an alternative embodiment of the present 
invention; and 



Figure 6 shows the frequency response of a filter with rectangular window of 
size N= 276, centered at 852 Hz, showing improved suppression of adjacent DTMF 
signal frequencies, in accordance with an alternative embodiment of the present 
invention. 



Figure 4 shows filter roll-off for a filter with rectangular window; 
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Figure 1 shows a typical tone detector. Pre- filter 1 filters the input signal to 
suppress out-of-band frequencies. The envelope detector 4 determines whether a 
signal is present or not and synchronizes the signal edges of the signal with the 
operation of the analysis filter bank 3. The analysis filter bank consists of multiple 
5 filters, one at each MF frequency of interest. These filters calculate the energy of each 
frequency. The energy calculator 5 computes the total energy of the signal after pre- 
filtering. The total energy combined with the results of the analysis filter bank 3 is fed 
to a decision logic block 7 that determines whether a valid tone is present or not. The 
decision block 7 may perform various tests, Such as minimum tone level, twist (if it is 
10 a dual tone), maximum frequency deviation, etc. As discussed in greater detail below, 
according to one aspect of the invention a new method is implemented in decision 
logic block 7 to determine frequency deviation from a nominal value by using the roll- 
off in the filter bank 3 combined with background energy. 

15 The analysis filters 3 are implemented as DFT filters using a modified 

Goertzel algorithm, one for each MF (e.g. R2) frequency. The Goertzel filter is 
centered at the R2 frequency. 

According to the preferred embodiment a rectangular window was chosen for 
20 the DFT (Goertzel) filter. The choice of a rectangular window results in reduced 
complexity relative to other filter types. For example, pre-processing the data with 
either a raised cosine or Kaiser window has a considerable MIPS impact. As 
indicated above, however, the leakage of the energies in the side lobes can be 
significant in rectangular windows (especially if two valid frequencies to be detected 
25 have a relative twist). Thus, according to an aspect of the invention, the filter size is 
chosen in such a way that the other MF frequencies fall in its spectral nulls, thereby 
greatly reducing the leakage effect. 
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R2 is a special case of MF communication. The six R2 signaling frequencies 
are situated at 540Hz, 660Hz, 780Hz, 900Hz, 1020Hz and 1 140Hz for the backward 
direction and at 1380Hz, 1500Hz, 1620Hz, 1740Hz, 1860Hz and 1980Hz for the 
forward direction. Thus, in each direction the frequencies are 120Hz apart. A window 
5 of size N has nulls situated at frequency intervals Fs/N, as shown in Figure 2. 
Consequently, for R2 signaling spectral nulls are needed at 120Hz intervals. 




In order to position the nulls at 120 Hz intervals, and assuming a sampling 
frequency Fs = 8kHz, a window size of 200 samples is chosen for the filter (the 
10 minimum window size for nulls at 120Hz is N = Fs / 120Hz, which is 66.67, such that 
the smallest integer window is thus N = 200). 



For example: If the frequency of interest is at 900Hz, the filter is centered at 
900Hz using a window of 200 samples. Spectral nulls are therefore located at each of 
15 the other MF frequencies, i.e. 540Hz, 660Hz, 780Hz, 1020Hz and 1 140Hz. Figure 3 
shows the R2 example for a frequency of 900Hz. It will be noted that the other R2 
frequencies fall exactly on the spectral nulls. 

For non-equally spaced frequencies (such as DTMF), the approach is not as 
20 straight forward as for the equally spaced frequency case. Nonetheless, the windows 
for each frequency detector (i.e. filter) may be sized in such a way that the closest 
frequency falls into a null. On the other hand, the window may be chosen such that the 
interference of the other frequencies is minimized. 
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For example, if the frequency of interest is the DTMF frequency at 852 Hz., 
the nearest adjacent frequency is 770 Hz (i.e. separated by 82 Hz). Choosing a 
window size of N=195 places nulls at intervals of 41Hz, as shown in Figure 5. The 
other DTMF frequencies (697 Hz, 941 Hz) do not fall on the spectral nulls, but are 
5 only suppressed by 24.8 dB and 23 dB respectively. However, if the window size is 
chosen to be N=276 as shown in Figure 6, none of the DTMF frequencies fall on 
spectral nulls. Nonetheless, the frequencies are suppressed by 25.3 dB, 24.8 dB and 
32.6 dB respectively, thereby improving the overall suppression. 



10 The filters comprising analysis bank 3 may be implemented in many ways. As 

indicated above, a Goertzel DFT filter was implemented according to the preferred 
embodiment. Implementations of Goertzel filters are well known in the art (e.g. 
Burrus, C.S. and T.W. Parks. 1985. DFT/FFT and Convolution Algorithms. New 
York: Prentice-Hall.). An example of a Goertzel filter for use in the multi-frequency 

15 tone detector of the present invention, centered at 900 Hz, is as follows: 

Fs = 8000; %Sampling Frequency 
Fc = 900; % Centre frequency 

coeff = (2*cos (2*pi* (Fc/Fs) ) ) ; % Goertzel coefficient 
N = 200; . % To create nulls at 4 0Hz intervals 



BEGIN 

% Feedback phase 
for k=l:N 

temp = input (k) + coeff *y(l) - y(2); 
25 y(2) = y(l) ; 

y(l) = temp; 

end 

% Feedforward phase (Calculates energy in Window) 
30 Yll = y(l) ; 

Y12 = y (2) ; 

energy = (Y11*Y11 + Y12*Y12 - Yll*Y12*coef f ) /N/N; 
return energy 

END 

35 

The return value "energy" is the tone energy at a specific frequency, which is 
fed into the decision logic 7. 
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The frequency deviation algorithm implemented within block 7 uses the roll- 
off of the filter 3 combined with an estimate of the background energy. Figure 4 
shows the roll-off characteristics of a rectangular window (for a size of 200 samples 
5 and a sampling rate of F s = 8 kHz). Thus, for a single tone frequency deviation of 1 5 
Hz from nominal, the roll-off is about -2 dB. In other words, if a tone deviates 15 Hz 
from its nominal frequency, its energy in the bin would be 2 dB down from its actual 
energy. This energy is not "lost", but has leaked outside the bin, contributing to the 
background energy. Consequently, by calculating the background energy the. 
10 frequency deviation of the tone may be determined. 

Total Energy is the energy of the output of the pre-filter 1. It is made up of the 
single frequency energy (E) and background noise (n). 

Total Energy = E + n 

15 The frequency energy is made up of the energy within the analysis filter 

(E_bin) and the leakage energy due to roll-off (E_ro). Thus, for a single tone: 

E = E_bin + E_ro 

If the frequency deviation is zero, then 

E = E_bin and Ero is zero 

20 Thus the Total Energy is made up of: 

Total Energy = E_bin + E_ro + n 

= E_bin + E_back 

where E_back is the background energy relative to the output of analysis filter 3. 
Thus, the background energy is calculated by subtracting the energy within the 
25 analysis filter from the total energy. The Signal-to-Noise Ratio (SNR) of the analysis 
filter energy relative to the background energy can be represented as: 

SNR ro = E bin/E back 
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= E_bin / (E_ro + n) 

Thus, there is a direct relationship between SNR__ro and E_ro, the leakage 
energy due to the roll-off, and consequently the frequency deviation. 

For a rectangular window: 

5 Roll-off = (sin(x) / x) 2 with 

x = 7i * f * size / Fs 

For tolerable frequency deviations up to f = 13 Hz, size = 200 and Fs = 8 kHz 
x = 1.021 and the roll-off is 0.6974. Assuming that the noise (n) is negligible relative 
to the leaked tone energy E_ro, then: 

10 E_back = E_ro = (1 - Roll-off)*E 

And E_bin = Roll-off * E 

SNR_ro = Roll-off / (1 - Roll-off) 

= 0.6974/(1-0.6974) 

= 2.3 

15 

SNR_ro provides an upper threshold limit for tolerable frequency deviation. If 
the ratio between E_bin and E_back is less than SNR_ro, frequency deviations of 
more than 13 Hz are accepted. If the threshold is set above SNR_ro, then frequency 
deviations of up to 13 Hz will not be tolerated. 

20 The dual tone case is very similar to the single tone case set forth above. 

However, in the case of dual tones the Total Energy is made up of both energies. 

Total Energy = El_bin + El_ro + E2_bin + E2_ro + n 

= El_bin + E2_bin + Eback 

And E_back = El_ro + E2_ro + n 

25 SNRl_ro = El_bin / (El_ro + E2_ro + n) 



• # 

-8- 



Assuming no twist between the tones (El = E2) and negligible noise, then for 
a frequency deviation of 13 Hz for both tones 

SNRl_ro = El_bin / E_back 

= Roll-off / ((1 - Roll-off) + (1 - Roll-off)) 

5 =1.15 

The situation is a bit more complicated if the tones have relative twist. Assume 

E2 = twist * El and twist = 6.3096 ( = 8 dB) with f = 13 Hz 

Then 

SNRl_ro = El_bin / E_back 
10 = Roll-off / ((1 - Roll-off) + twist * (1 - Roll-off)) 

= Roll-off / ((1 + twist) * (1 - Roll-off)) 
= 0.315 
= -5dB 

In this case, SNR_ro provides the upper threshold limit which, it will be noted, 
15 is considerably relaxed relative to the threshold without twist, in order tolerate a 
frequency deviation of at least 13 Hz when dual tones and twist are present. 

The following is a pseudo code implementation for a typical dual tone detector 
using roll-off: 

BEGIN 

20 Calculate Total Energy 

Calculate bin energies at each nominal frequency- 
Select two highest frequency bins, ie . El_bin and E2_bin 
Calculate E_back = Total_Energy - El_ bin - E2_bin 
If El__bin < SNRl_ro * E_back 

25 Status = failed 

Elseif E2_bin < SNR2_ro * E_back 
Status = failed 

End 

END 

30 



Although a preferred embodiment of the present invention has been described, 
those skilled in the art will appreciate that variations and modifications may be made. 
For example, the embodiment set forth above is implemented via Goertzel DFT 
filters, whereas the filters 3 can also be implemented as FIR filters or IIR filters. This 
and all other such alternatives and variations are believed to be within the sphere and 
scope of the invention as defined by the claims appended hereto. 



